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We demonstrate that the instabilities of asymmetric nuclear matter at sub-saturation 
densities do not present two types of instabilities as usually discussed but a unique one. 
The associated order parameter is everywhere dominated by the isoscalar density and so 
the transition is of liquid-gas type even in the so-called chemical instability region. 
| Below saturation density, the nuclear interaction is expected to lead to a liquid-gas 
phase transition [ P. Recently, a converging ensemble of experimental signals seems to 
have established the phase transition. One is the spinodal decomposition [|2j. It has been 
argued that asymmetric nuclear matter do not only present a mechanical instability for 
which the density is the order parameter but also a broader chemical instability associated 
with fluctuations of the matter isospin content [E]. Indeed, it is usually argued that it 
exists a region in which the compressibility at constant asymmetry and constant isospin 
is negative (see figure (JIJ) leading to the interpretation that the system is mechanically 
unstable. Above a maximum instability the isotherms at constant asymmetry does not 
presents any back bending leading to the idea that the system is mechanically stable. 
However, looking at the equilibrium of the chemical potentials one can see that above this 
maximum asymmetry for mechanical instabilities the system may amplify fluctuations in 
the proton neutron concentration leading to a second instability region usually called 
chemical instabilities. 

However, we have recently shown that this splitting of the spinodal region in to two 
types of instabilities a mechanical and a chemical one is not correct and that ANM present 
only one type of instabilities [ |3] • The associated order parameter is dominated by the 
isoscalar density and so the transition is of liquid-gas type. The instability goes in the 
direction of a restoration of the isospin symmetry leading to a fractionation phenomenon. 
Our conclusions are model independent since they can be related to the general form of 
the asymmetry energy. They are illustrated using density functional approaches derived 
from Skyrme and Gogny effective interactions. 

Let us consider ANM characterized by a proton and a neutron densities p, = p p , p n . 
These densities can be transformed in a set of 2 mutually commuting charges p% — p±, p% 
where p\ is the density of baryons, p\ = p n +p p , and p3 the asymmetry density ps = p n —p P - 
In infinite matter, the extensivity of the free energy implies that it can be reduced to a 
free energy density : F(T,V, Ni) = VJ r (T,p i ). The system is stable against separation 
into two phases if the free energy of a single phase is lower than the free energy in all 
two-phases configurations. This stability criterium implies that the free energy density is 
a convex function of the densities pi. A local necessary condition is the positivity of the 
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Figure 1. The pressure as a function of the total density for a fixed temperature and 
different asymmetry. The region where the isotherms present a back bending are usually 
related to a mechanical instability. However, one can define a broader region of instability 
by allowing the isospin fluctuation (traditionally discussed as a chemical instability). We 
will show that this classification is artificial since it does not correspond to a modification 
of the physical properties of the instabilities. 

We show in the left part of Fig. El the spinodal contour in ANM for several forces. 
It exhibits important differences. In the case of SLy230a force (as well as SGII, DIP), 
the total density at which spinodal instability appears decreases when the asymmetry 
increases whereas for SIII (as well as Dl, D1S) it increases up to large asymmetry and 
finally decreases. We observe that all forces which fulfill the global requirement that they 
reproduce symmetric nuclear matter (SNM) equation of state as well as the pure neutron 
matter calculations, leads to the same curvature of the spinodal region. We can appreciate 
the reduction of the instability when we go away from SNM. However, large asymmetries 
are needed to induce a sizable effect. The temperature dependence of the spinodal contour 
can be appreciated in the right part of Fig. |21 As the temperature increases the spinodal 
region shrinks up to the critical temperature for which it is reduced to SNM critical point. 
However, up to a rather high temperature ( 5 MeV) the spinodal zone remains almost 
identical to the zero temperature one. 

In the considered two-fluids system, the [J 7 ^] is a 2 * 2 symmetric matrix, so it has 2 
real eigenvalues A ± [|7j : 
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Figure 2. Projection of the spinodal contour in the density plane : left, for Skyrme 
(SLy230a SGII [E], SIII and Gogny models (Dl [Q2|, D1S [Ej, DIP [d]) 
; right, temperature dependence of the spinodal zone computed for the SLy230a case. 



associated to eigenvectors dp defined by (i ^ j) 

Spf _ F ij _ ^ — ^ 
8pf A± - :F) , Tij 

Eigenvectors associated with negative eigenvalue indicate the direction of the instability. 
It defines a local order parameter since it is the direction along which the phase separation 
occurs. The eigen values A define sound velocities, c, by c 2 = j^Pi A. In the spinodal 
area, the eigen value A is negative, so the sound velocity, c, is purely imaginary and the 
instability time r is given by r = d/\c\ where d is a typical size of the density fluctuation. 

The requirement that the local curvature is positive is equivalent to the requirement 
that both the trace (Tr[.Fjj] = A + + A~) and the determinant (DetfjF^] = A + A~) are 
positive 

Tr^] > 0, and DetfJ^-] > (4) 

The use of the trace and the determinant which are two basis-independent characteristics 
of the curvature matrix clearly stresses the fact that the stability analysis should be 
independent of the arbitrary choice of the thermodynamical quantities used to label the 
state e.g. (p p , p n ) or (p u p 3 ). 

A different discussion can be found in the literature [EJ|3] and we will now clarify the re- 
lation of this discussion with the eigen modes analysis. Indeed, from the thermodynamical 
relation 

PnDetlFij] = ^\t,p J^|t,„ (5) 

one can be tempted to relate separately ^-\t,p and §^\r, y to the two eigenvalues A + and 
A~. The discussion of the sign of A + and A~ is replaced by the discussion of the sign of 
^j-\t,p (also called chemical instability) and §^\r,y (also called mechanical instability). 
This is correct in SNM and the mechanical instability is associated with density fluctu- 
ations whereas the chemical instability is associated with the isospin density fluctuation. 
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In the case of ANM, because natural symmetries are lost, theses relations break down 
where one can show that [|7] 

dP A + A~ 

— — \ T>y = — F (t cos/5 + sin/5) 2 H — j=(t sin/3 — cos/5) 2 (6) 
opi ' \/t yt 

1 

dp p 



dy 




p n A + A~ - — \ T ,y (7) 



where /3 = 1/2 tan J- np j{T pp — J- nn ) and * = PhNq / p p N$. These two equalities illustrate 
explicitly that the simplicity of SNM is not preserved in ANM. Hence, one should come 
back to eigen analysis of the curvature matrix [J 7 ^]. 

The local stability condition in ANM is express by Eq. HJ If is violated the system is 
in the unstable region of a phase transition. Two cases are then possible : i) only one 
eigenvalue is negative and one order parameter is sufficient to describe the transition or ii) 
both eigenvalues are negative and two independent order parameters should be considered 
meaning that more than two phases can coexist. 

For ANM below saturation density, the case ii) never occurs since the asymmetry energy 
has always positive curvature (^33) ■ Indeed, the asymmetry term in the mass formula 
behave like (N — Z) 2 times a positive function of A showing that the dominant p% depen- 
dence of the asymmetry potential energy is essentially quadratic and that JF 33 is a positive 
function of the total density. Recent Bruckner calculations in ANM [ [[| have confirmed 
the positivity of JF 33 . They have parameterized the potential energy with the simple form 
V(pi,P3) — Vo(pi)pi + Vi(pi)p| with V0/V1 ~ —3 and Vo < 0. This is also true for 
effective forces such as Skyrme forces. For example, the simplest interaction with a con- 
stant attraction, to, and a repulsive part, t 3 pi, leads to V(pi,p 3 ) = (3p 2 — p 2 )/5(pi) with 
/5(pi) = (to + ^3/6pi)/8. The function /3(pi) is negative below saturation density, hence 
the contribution of the interaction to .F33 in the low density region is always positive. 

These arguments show that, below saturation density, the P3 curvature, JF 33 , is expected 
to be positive for all asymmetries. Since the curvature in any direction, J-'a, should be 
between the two eigenvalues A~ < Ta < A + we immediately see that if JF 33 is positive one 
eigen curvature at least should remain positive. In fact for all models we have studied J-33 
appears to be always large enough so that the trace is always positive demonstrating that 
A + > 0. Since Trfj^V,-] = T nn + J~ pp , this can be related to the positivity of the Landau 
parameter T nn and T m . 

The large positive value of J-33 also indicates that the instability should remain far 
from the p 3 direction i.e. it should involve total density variation and indeed we will see 
that in all models and for all asymmetry the instability direction hardly deviates from a 
constant asymmetry direction (<5p 3 <C dpi). This isoscalar nature of the instability can be 
understood by looking at the expression of the eigen-modes in the (p„; p p ) coordinates. 
Since A~ < Ta < A + , the differences A~ — Ta is always negative demonstrating, using 
Eq. El that the instability is of isoscalar type, if the Landau parameter T np is negative. 
Hence, there is a close link between the isoscalar nature of the instability and the attraction 
of the proton-neutron interaction [|7j. 

Various contour of equal imaginary sound velocity are represented in Fig.Elfor SLy230b 
and DIP interactions. The more internal curve correspond to the sound velocity i0.09c, 
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Figure 3. This is the projection of the iso-eigen values on the density plane for Slya 
(left) and DIP (right). The arrows indicate the direction of instability. The mechanical 
instability is also indicated (dotted line). 

after comes z0.06c, i0.03c and finally 0, the spinodal boarder. For these two recent forces 
which takes into account pure neutron matter constraints, the predicted instability do- 
mains are rather similar. We observe that in almost all the spinodal region the sound 
velocity is larger than 0.06c. 

Let us now focus on the direction of the instability. If 5p~ is along y=cst then the insta- 
bility does not change the proton fraction. For symmetry reasons pure isoscalar (<5p 3 = 0) 
and isovector (5pi = 0) modes appears only for SNM so it is interesting to introduce a 
generalization of isoscalar-like and isovector-like modes by considering if the protons and 
neutrons move in phase (5p~5pp > 0) or out of phase (5p~5pp < 0). Fig. El shows the 
direction of instabilities along the spinodal boarder and some iso-inst ability lines. We 
observed that instability is always almost along the pi axis meaning that it is dominated 
by total density fluctuations even for large asymmetries. In fact, the instability direction 
is always between the y=cst line and the pi direction. This shows that the unstable direc- 
tion is dominated by isoscalar modes as expected from the attractive interaction between 
proton-neutron. The total density is therefore the dominant contribution to the order 
parameter showing that the transition is between two phases having different densities 
(i.e. liquid-gas phase transition). The angle with the pi axis is almost constant along a 
constant y line. This means that as the matter enters in the spinodal zone and then dives 
into it, there are no dramatic change in the instability direction which remains essentially 
a density fluctuation. Moreover, the unstable eigenvector drives the dense phase (i.e. the 
liquid) towards a more symmetric point in the density plane. By particle conservation, 
the gas phase will be more asymmetric leading to the fractionation phenomenon. Those 
results are in agreement with recent calculation for ANM [|7j and nuclei [ITU]. 

We want to stress that those qualitative conclusions are very robust and have been 
reached for all the Skyrme and Gogny forces we have tested (SGII, SkM*, RATP, Dl, 
D1S, DIP,...) including the most recent one (SLy230a, DIP) as well as the original one 
(like Sin, Dl). 

In this paper, we have shown that ANM does not present two types of spinodal insta- 
bilities, a mechanical and chemical, but only one which is dominantly of isoscalar nature 
as a consequence of the negativity of the Landau parameter T np . This general property 
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can be linked to the positivity of the symmetry energy curvature JF 33 . This means that 
the instability is always dominated by density fluctuations and so can be interpreted as a 
liquid-gas separation. The instabilities tend to restore the isospin symmetry for the dense 
phase (liquid) leading to the fractionation of ANM. We have shown that changing the 
asymmetry up to p p < 3p n does not change quantitatively the density at which instability 
appears, neither the imaginary sound velocity compared to those obtained in SNM. All 
the above results are not qualitatively modified by the temperature which mainly intro- 
duce a reduction of the spinodal region up to the SNM critical point where it vanishes. 
The quantitative predictions concerning the shape of the spinodal zone as well as the 
instability times depends upon the chosen interaction but converge for the various forces 
already constrained to reproduce the pure neutron matter calculation. 
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